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A CLASS OF PERIODIC SOLUTIONS OF THE PROBLEM OF THREE 
BODIES WITH APPLICATION TO THE LUNAR THEORY* 

BY 

FOEEST EAY MOULTON. 

1. Introduction. 

2. The Differential Equations. 

3. Existence of the Periodic Solutions. 

4. Practical Construction of the Solutions. 

5. Applications to the Lunar Theory. 

6. Applications to Darwin's Periodic Orbits. 

1. Introduction. 

This paper treats certain periodic solutions of the problem of three finite 
masses which move in the same plane under their mutual attractions. The 
forces are supposed to vary according to the Newtonian law, and it is assumed 
that the bodies are spheres. The periods of the solutions are equal to the 
synodical periods of the bodies, and the orbits reduce to circles when the dis- 
turbing forces are put equal to zero. It follows that they are of the class 
called by PoiNCARjfi " Solutions cle la premiere sorte." -f 

In the series which have been developed for expressing the coordinates of the 
moon there are certain terms whose period is the same as that of its synodical 
revolution. In the notation of DelaunayJ they are the terms which involve 
sin ^D and cos JcD , where k is an integer. It is shown that these terms belong 
to the periodic orbits treated in this paper. By putting the mass of the moon 
equal to zero and transforming the parameter in terms of which the expansions 
are made, we can compare the coefficients of the various terms with those given 
by Delaunat.§ 

As applied to the Lunar Theory the results obtained in this paper are, when 
the mass of the moon is neglected, coextensive with those given by Hill in his 
Iiesearches,\\ and by Brown in his first memoir.^ 

* Presented to the Society under a different title, September 8, 1905 and April 14, 1906. 
Eeceived for publication June 26, 1906. 

t Les Metliodes Noiwelles de la Mccanique Celeste, vol. 1, p. 97. 
t Memoires de I'Aoade mie des Sciences, vol. 29 (1867), p. 801. 
§ Loo. cit., pp. 803-924. 

II American Journal of Mathematics, vol. 1 (1878), pp. 5-26, 129-147, 245-260 
If American Journal of Mathematics, vol. 14 (1892), pp. 141-160. 
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538 F. E. MOULTON : THE PROBLEM OP THREE BODIES [October 

When the masses of two of the bodies have the ratio 1 to 10, while that of 
the third is infinitesimal, the problem of this paper reduces to that which 
Darwin treated by numerical processes.* For the sake of comparison, one of 
his " Satellite A " orbits and one of his " Planet A " orbits have been computed 
from the literal series derived below. A retrograde satellite orbit having the 
same sidereal period is also given. Darwin's " Satellites B " and " Satellites 
O" are imaginary for small values of the disturbing forces. They belong to 
values of the parameter which lie beyond that for which the series have a singu- 
larity, and their coordinates, therefore, are not represented by these series. 
There are only three classes of orbits of the type in question, real or imaginary, 
in which the bodies revolve in the forward dii*ection, and it appears, therefore, 
that Darwin's search for these periodic satellite orbits was exhaustive. How- 
ever, the number of retrograde orbits is the same as the number of those which 
revolve in the forward direction. 

After having selected variables which seem most convenient for the construc- 
tion of the series and for practical use, the existence of the solutions of the type 
in question is established. For this purpose the process of analytical continua- 
tion is employed which Poincar6 developed in this general connection, and 
which he made the central idea of his researches on the figures of equilibrium 
of rotating fluid masses f and on the problem of three bodies.^ After the 
demonstrations of the existence and the properties of the solutions have been 
made, the practical construction of the series is carried out along lines which 
have been shown to be legitimate. Instead of arranging the solutions as Fourier 
series, as has been customary among astronomers, power series are employed, 
and the work has the simplicity which is characteristic of processes involving 
power series. For example, every coefficient is determined by a single step and 
is modified by no subsequent operations. Moreover, it is possible to prove the 
convergency of the series within certain limits, though the precise limits where 
convergency ceases have not been found. 

2. The differential equations. 

Let the three masses be £J, M, and S. Choosing ^ as the origin and 
letting the rectangular coordinates of JTand S he x, y and ^, Irrespectively, 
we have the differential equations of motion : 



*Aota Mathematica, vol. 21 (1897), pp. 99-242; Mathematische Annalen, vol. 51 
(1898-9), pp. 523-583. 

t Acta Mathematioa, vol. 7 (1885), pp. 259-380. 

J Bulletin Astronomiqae, vol. 1 (1884), pp. 65-74; Acta Mathematioa, vol. 13 
(1890), pp. 8-270 ; Les MSfhodes Nouvellet de la Mecanique Celeste (1892-9). 
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where 



rPV V I V— ?/ 7/ 



T = Vx^ + y\ p=V{x-Xf + {y-rf, R = vx^ + r\ 

We shall now refer the motion of S to the center of gravity of E and M by 
the transformation 

Then passing to polar coordinates by the transformation 

X ^ r cos « , x' ^ r' cos «', 

2/ = rsin«, 2/'^'*'*^'*''''' 

the differential equations become 



df ^ 



(S)V'^^^^=^=«(^-^)»*- 



(1) 



'^^^^dt-dt=^^^^^^yr^§.^ 

d'^r ^,(dv\ ¥E{E+M+S) _¥{E + M+ S)dv^ 



\dt) 



where 



df \dt)^ (E+My^ E + M dr 

, dH' dr' dv _ ¥{E+M+ S) 1 du^ 
^ ^df"^ 'di di~ E~^M / W ' 






Trans. Am. Math. Soe. 36 
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If rjr < 1 then IjH — l/r and 1/p caa be expanded as converging power 
series in r/r for all values of v and v. Only cases in which this inequality is 
satisfied will be treated. In the Lunar Theory r/r is about ^^^, and in the 
satellite orbits computed by Darwin, about Jj. 

The series for u and u are 

1 f 1 1 B + M /rW _ ^ 

(E'^ — 3P)fr V 
+ {WT'Mf Vr' j * t^ ""^ ^^ - «') + ^ «°' 3 (^■ - ^ )] 

+ r^-+-Jfy(.')6V[9 + 20cos2(^-«') + 35cos4(.-^')] 
+ r^ijrMY ( 7 ) TbP^«oK'"-'"')+35<5o^^(«--''')+63cos5(«-v')] 



+ 



« =7P^+-~4"(7J [l + 3cos2(t-«')] +.. 
With these values of « and ii' equations (1) become 

+ f [-J-^Jp [3 cos (,;_«') + 5 cos 3(^ - ^' )] + ■•• }, 
cPv ^dr dv Je^JS r \ „ . „ , ,. 



f ^ — M) r 
(2) + I |^^:p-^^ 7 [sin (t, - t,' ) + 5 sin 3 (t, - t;' )] + 



iV ,fdv'y ¥{E+M+8)_ ¥EM{E + M+_S) 



_ ,fdv'Y _ 

5<^ ~ '■ \dt ) ^ 7'' ~ ~ 2(E+MY 

r^ I 
X -,, [l + 3cos2(«-^'')] + 



,d'v' „dr' dv' PEM{E+M+ S) ^ , „ . „/ 
'"df+^ Tt dt = "2X^+1/7 ,- i 3 sm 2 (. - « ) + 
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When the right members of these differential equations are neglected, that is, 
when M is supposed to revolve around E without having its motion disturbed 
by S and when S is supposed to revolve around the center of gravity of E and 
J/^ without having its motion disturbed because of the distribution of the mass 
between the two separate bodies E and M, a particular solution of the equa- 
tions is 



(3) 



hVE + M 
V = 5 t = nt. 



, kVE + M+S 
r ■= a , w = 5 i = n t. 



where n and n' are the angular motions of E and S respectively. In develop- 
ing the periodic solutions from the differential equations including the right 
members, the mean angular motions will be kept the same as in the undisturbed 
orbits, and n and n will be regarded as constants which are given by observa- 
tion (in the applications to the Lunar Theory), or which are assumed. The 
masses are regarded as given, and the constants a and a are defined by the 
equations 

(4) n'a?=¥{E+M), 7i\'' =¥{E + M + S). 

New variables^, w. p' and w' will be defined by the equations 

r = a(^l + p), V = nt + w, 

(5) ,,,,,, 
r =a'(l+p'), v' = n't + w'. 

It follows that ap,w, a'p and w are the deviations from the circular orbits of 
M and 8 due to their mutual perturbations. To reduce the differential equa- 
tions to a still more convenient form we shall let 

n n {\ -\- m\ 
-"' 7 , whence n = — ^ , 



n — n 



71 

(6) r=(n — n')t=— t, 

^ ' ^ ' m 

S^ 

'^- E+M+ S- 

It foUov/s from (4) and the definition of m that 

E+ M 



a / m y/ 
^ ' « ' ^ \ 1 + m ) \ 



E+M+ 8 J 



That is, the ratio aja is expressible in terms of the parameter m and might be 
eliminated from the equations. But in the Lunar Theory aja' is regarded as 
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a separate parameter and is retained explicitly. In order to keep the results in 
such a form that comparison may be made with the work of Dblaunay, and 
at the same time to avoid the inconvenience of having power series in two 
parameters, we shall put in general 

^ ^ {E + My \a J -'-1 ' \a ) ' ' 

where the M. and iV^. are determined so that (7) shall be fulfilled. The m will 
be regarded as a variable parameter in the argument, and the M. and N. as 
fixed constants. Therefore the m in (7) and the m in (8) belong to different 
problems except when they give the same value of a\a' . At the end of the 
work this relation will always be fulfilled, and the equations (8) amount merely 
to an allowable generalization of the meaning of m . The relations have been 
selected so that the 3f. and If. are approximately (f )■' in the case of the 
Lunar Theory, and the power of m indicates the order of the term. With these 
substitutions the differential equations become 

1 + p 
+ |J/jm^ , [3cos(t + w — w') + 5 cos S (r + w — w')] 

(1 +pf 
+ IM^m" ":^ 'Y [9 + 20 cos 2{t+w-w') + S5 cos 4{t+w-w')] 

(1 + nf 
+ gSjJIfjm^pj— J-^3 [30cos(t + w — w')4- 35 cos 3(t+ w — w') 



+ 63 cos 5(t + w — w')] + • • 

1 + » 
(9) + |J/j7n^3 -, [sin(T + w) — w') + 5 sin3(T4- w — w')] 

+ |J!f^m^ ^ ^; 2 [ 2 sin 2 ( T + w — w') + 1 sin 4(t + m — w')] 

(1+pf 
+ -i|Jf,m^^-j ^, [2 sin(T + w — w') + 1 sin3(T + w — w') 



+ 21 sin 6 ( T + w — t«') ] -I i , 



1906] WITH APPLICATION TO THE LUNAR THEOEY 543 

X I [1 + 3cos2(t+ w — w')] + 
,. ,^(Pw' ^ dp' f dw' \ 



^^^._^^eil±Pl 



2(U+M)-'' {1+p'y 

X ]3sin2(T + w — w') + 



These equations are general for the determination of the motion of three 
mutually attracting spheres except for the condition that a ( 1 + p) shall be less 
than a'(l + p') throughout the motion. 

3. Existence of the periodic solutions. 

When m = equations (9) admit the periodic solution p = w=^' = w'=0. 
It will now be proved that for ?» 4= 0, but sufficiently small, the equations 
admit a periodic solution expansible as power series in m and vanishing with ni . 
It is the analytic continuation of the solution p =^ w= p' ^ w' = m ^ with 
respect to m as the parameter. Since t enters explicitly in (9) under the sine 
and cosine functions in integral multiples, it follows that the period of the solu- 
tion must be a multiple of 27r . 

It can easily be shown from the vpay in which p, w,p', w', and t enter in 
(9) that if we impose the conditions dp/dr = w = dp' /dr = w = at t =^V, 
then p and p' are even functions of t — j-tt, and w and w' are odd functions of 
T — JTT. This means geometrically that at t = jV the bodies are in a line and 
are crossing it perpendicularly, and that the orbits of M and S are respectively 
symmetrical with respect to axes coinciding with this line at r=JTr and rotating 
around E and the center of gravity of E and M at any constant rate. Orbits 
of this type, whether they are periodic or not, will be called symmetrical orbits. 

The existence of symmetrical periodic orbits will be established first. Sup- 
pose that, at T = , 



r = a(l +|j) = a(l-f a)(l — e), /= a'{l -\- p) = a{l + a')(l — e), 

dr dp dr ,dp' 

dr ~ dr ~ '' dr ~ dr ~ '' 



(10) 

W = , M)' = , 



dw (l+m)vl + e dio mV^ + e 
1 + "* + 77- = > V- ,— -\l7-i rva' "^ + ^jZ= 



dr -(l-j-a)3(l-e)S' ^ dr " (1 -j- «')?(!_ e')?' 

If the right members of equations (9) are neglected, and if we put ?• = a (1 +|)) , 
( 1 -f- m ) T + M) = r< , r'= a'( 1 + 1)' ) , mr -f w'= ii, then the left members reduce 
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to two sets of equations for the two body problem. It is found that, with the 
initial conditions (10), aa, and a a are the increments to the semi-axes of the 
orbits of M and 8 respectively, and that e and e are the eccentricities of their 
orbits. Because of the well-known properties of the solutions in terms of these 
parameters, the properties of the general solutions so far as they do not depend 
upon the right members of (9) can easily be written down. For this reason 
these particular parameters were selected. 

Equations (9) are to be integrated as power series in a, e, a', e and m, van- 
ishing with these parameters. It follows from their form and Poincae^'s 
extension * of Cauchy's theorem that, when an arbitrary interval for t has been 
chosen in advance, these parameters can be assigned such smaU moduli that the 
series will converge uniformly as t takes all values in the interval. We shall 
choose as the interval • • • 2^7r (k an integer) and the solution may be written 

p = |)j(a, 6, a', e', to; t 
dp 



di 



= p^{a, e, a , e ,ni; t 



(11) 



w = p^(^a, 6, a', e', m; r 

dw ^ , , 

-^ =^^(a, e, a, e , m; t 

^' = |)5(a, e, a', e', m; r 



dw' 

dT 



Pg(a, e, a', e', w; t), 



where p^, • • • , !f3g are power series in a, e,a', e' , and m . The independent vari- 
able T enters in the coefficients. 

It follows from the symmetry of the orbits under the initial conditions (10) 
that if at T = kir the three bodies are again in a line and are crossing the line 
perpendicularly, then the orbits necessarily reenter at t = 2^77, and conversely. 
Hence the necessary and sufficient conditions for a symmetrical periodic solu- 
tion of period 2kTr are 
(12) p.(a, e, a', e', m; ^tt) = (j = 2,3,6,7). 

Since ^2 = ^3 — ^6 = ^7 = at t = , it follows that equations (12) are iden- 
tically satisfied hja=e = a' = e'=m = 0. It will be shown that they can 
be solved for a, e, a' and e' as power series in m, vanishing with ni , and con- 
verging for sufficiently small values of m . When the results of these solutions 



*Acta Mathematioa, vol. 13 (1890), p. 16; Les Methodes Nouvelles de la Mecanique 
Celeste, vol. 1, p. 58. 
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are substituted in (11), the coordinates will be expressed as power series in m 
alone, and they will converge for sufficiently small values of m . The conditions 
for periodicity will have been satisfied, and therefore for m sufficiently small, 
but not zero, the series will converge uniformly for all values of t. Since the. 
series converge for a continuous range of values of m , the coefficie7it of each 
power will he periodic. 

The first two equations of (9) do not carry p' and w' except in their right 
members ; therefore ^)j, • • •, p^ (io not contain a and e' except in terms which 
are multiplied by m^. Likewise p^, ■ ■-, pgdo not involve a and e except in 
terms which are multiplied hy m^ . All the tei'ms in pj , • • • , |3^ which do not 
contain m^ as a factor can be obtained from the- left members of the first two 
equations of (9) set equal to zero. All the terms in p^, • • •, p^ which do not 
contain m^ as a factor can be obtained from the left members of the third and 
fourth equations of (9) set equal to zero. The left members of (9) set equal 
to zero constitute two pairs of equations of the two body problem, and making 
use of this fact and the initial conditions (10), we may write the solution (11): 

2> = pj= — e {cos i'T + 2(cos2i'T — 1) + . • •} -\-m^'p[{a, e, a', e', m; t); 

3^ = p2 = e/^lsin VT + esin2i/T + .. .} -f ni'')(:>'^(a, e, a', e', m; t), 
w=1()^= { — (l+m)T + vT+2esinw+---} + m^'p'^{a, e, a', e',m; t) 
-^ = Pi= {—(t + m) + V -i-2eveosvT+ ■■■} + m%(a, e, a', e', m ; t), 
p' = p^=z — e'{ cos v't + 2" (cos 2i/'t — 1) H } + mY^{a, e, a, e', m; t), 

-|--= pj. == e'v'{sin v't + e' sin 2k't + • ••} + m^p'^^a, e, a', e\ m; t), 
w' =pj= { — mr -{■ v't H- 2e' sin i/'t -f- • ••} + m^\)'^{a, e, a', e', m; t), 

dw' , , n I , , ■. c , , , , 

-J- = ^8= { — ^ + I* + 2e 1/ cosv T + • . •} + TO''pg(a, e, a , e , m; t), 

where 

(14) .= ! + '" 



(! + «)«' -(i4.«')i 

The terms contained in pj , • • ■ , p^ depend upon the right members of the 
differential equations. 



(15) 
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The conditions for periodic solutions (12) become 

^^^(^Tr) = et^jsin vh-rr + e sin 2vhir +•••}+ mh;)'.J^a, e, a , e', m ; Icir) = , 
^3(^77) = { — (1 + m)h-K + vh:: + 2e sin vkir + • • •} 

+ m^Tf'J^a, e, a , e', m; kir) = 0, 

'^Jyh-jr) = e'i'fsin i''i7r+e' sin ^v'h-7r-\ } + »i''pg(a, e, a', e', m; Ajtt) = 0, 

J3j,(/<;7r) = { — mk-K + v'kir + 2d?' sin i''^7r + • • •} 

+ m^Py(o(, e, a', e', m; ^tt) = 0. 

All the unwritten terms in the { | involve sines of multiples of vkir or v'kir ; 
those in the first two { } carry e^ as a factor, and in the last two, e' . 

Substituting the expression for v' as defined in (14) in the last two equations 
of (15), we find that ^^(^Tr) is divisible by m^ and that ^^(^-77) is divisible by 
m. After division by these factors there is no term independent of a\ e', and 
m, and the determinant of the linear terms in a and e' is 



A = 



0, k-JT \ 

\ 

— ^k-K, Ik-TT \ 



which is distinct from zero. An exception can occur only if mk is an integer, 
when the determinant does vanish. But since k may be put equal to unity this 
requires that m shall be an integer, a case which arises for no satellite and for 
which the series derived are almost certainly always divergent. Heuce it will 
be supposed that m is not an integer. Since A is not zero, pg(/i,-7r)/m^ = and 
}jj (i7r)/m^ = may be solved uniquely for a and e' as power series in a, e, and 
m. These equations are satisfied by a' = e' = lyi = whatever a and e may be, 
and therefore the solutions contain m as a factor, and may be written 

a' = m5',(a, e, m), e' = mq^[a, e, m). 

When these expressions for a and e are substituted in the first two equations 
of (15) the latter become 

b„ ( kir ) = ei/ 1 sin vk-rr + e sin Ivk-rr +•••} + wi^ p" ( a , e , m ; kir ) = , 
(16) 

^,(A;7r) = [~ (1 -\-m)k'7r+ vkir +2e sin vk-rr + • • ■} +-m^T(i"^{a, e, m; A;7r)=0. 

Before solving equations (16) we shall find the degree of the lowest terms in 
m alone which they contain. The terms which are independent of a, e, a' and 
e' are found by integrating equations (9) in series of the form 

OO « CC CC 
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Substituting these expressions in (9), equating coefficients of corresponding 
powers of m, integrating the resulting equations, and imposing the initial con- 
ditions (10), we find 

Pj = Wj = , 

p^ = — Tj -j- 27) cos r — rj cos 2t , 

w.^ = ^rjT — 4r) sin T + -^-7; sin 2t , 

j)^ = I j; 4- ^T] cos T — ^T) COS 2t — 2ijT sin t , 

Wj = 7;t — ^rj sin t + J^^-tj sin 2t — 4jjt cos t . 

These expressions give, at t = Ictt, 

dp. dj) 



^i'3 



fZ' 



= — (— 1 )" 2T)h-ir , Wj = 7?^7r — {—If 4:r)k'7r . 



Therefore the lowest power of m contained in the first equation of (16) is the 
third, and in the second equation, the second. 

Substitute in (16) the expression for v given in (14) and consider the solu- 
tion of the second equation. Every term contains either m or a as a factor. 
The lowest degree in a alone is the first, and in m alone, the second. Therefore 
the equation can be solved for a as a power series in m and e, containing m^ as 
a factor. Thus 

(17) a = rn^q^^e, m) . 

Let (17) be substituted in the first equation of (16). After the substitution 
it is found that a single factor m can be divided out leaving a linear term in e 
and a second degree term in m as the lowest. Consequently the equation can 
be solved for e as a power series in m , containing m^ as a factor. Therefore 
we have 

(18) e ^ m^ q^^m) . 

By substituting back through the several steps of the solution, we obtain 

a = m''q[{m), e = m^q'^{m), 

a ^= 7n? q'^{m) , e = m? q\{in) , 

where q'l, ■ ■• iq^ are power series in m . The first two contain terms independent 
of m ; the question for the last two will be left open until the practical construc- 
tion of the periodic series is considered, though it might be answered at this 
place. 
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When the expressions (19) are substituted in equations (13) the periodic solu- 
tion is obtained. The series Jbr p and w start with terms of the second degree 
in in. When k has been chosen the whole process is unique. Therefore for a 
given value of m there is but one direct orbit of the type under consideration 
which reenters after k synodic revolutions. This holds when k=l. But all 
the reentrant orbits after k revolutions include those which reenter after one 
revolution ; therefore, since the result is unique for every k , all of the orbits of 
the type treated reenter after a single synodic revolution. 

If the orbit of JIf is retrograde, n is negative and m has the definition 

— n 
m ^ 



n + n 

With this definition of m the discussion of the existence of symmetrical periodic 
orbits proceeds precisely as before, and it is found that they exist. Conse- 
quently, there is but one direct and but one retrograde orbit of the type in 
question with a given sidereal jjeriod. For the retrograde orbit the modulus 
of m is smaller than for the direct, and the deviations from a circular orbit are 
in general less. 

Suppose the condition is no longer imposed that the orbits shall be symmet- 
rical, and consider the question of the existence of periodic solutions. We may 
suppose that at t = 

dp dw 

p^a, w=0, ^=^, -^=y. 

That is, at the origin of time £J, M and S are in a line, but the initial condi- 
tions are not otherwise defined. The conditions w = w' = are no limitations 
upon the generality of the problem for they amount simply to a particular 
choice of the origin of time. The solutions may be expressed as power series in 
a, /3, 7, a', /S', 7' and m whose moduli may be taken so small that the conver- 
gence holds for T in the whole interval • • • 27r . The conditions for periodic 
solutions with period lit are 

p{a, ^, 7, a', /3', 7'; 27r) =;?(«, ^, 7, a', ^', 7'; 0) = a. 



(20) ^(a,/S,7,«',/3',7'; 2,r) = ^ (a, /3, 7, «', yS', 7'; 0) = /3, 



dw' , - / ^, / r, ^ ^''■'^' , r. , ^, , r. ■. , 

-^{a, ^,y,a,^,y; 27r) = -^(a, /3, 7, a, ^, 7; 0) = 7. 
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Equations (9) admit two uniform integrals, tliat of vis viva and that of areas. 
They may be written 

f dp dw , dp , dw'\ , n' I N 

„ / dp dw , dp , dw \ -r^ , n ' n' ' ^ 

^- Y^ dr^ ^' Jt'^^' ^ ' ^' 1!? j = ^^(^' ^' 'y' ^' ^''y' "')• 

It follows from these relations that when six of equations (20) are satisfied the 
remaining two are reduced to identities.* Consequently two of the equations 
may be suppressed and the remaining six solved for a, /3, 7, a, /3', 7' in terms 
of m. These solutions furnish the periodic orbits. 

In order to determine which equations may be most conveniently suppressed 
because of the existence of the integrals i^j and F^ some details are necessary. 
Likewise in solving the remaining six equations more suitable parameters must 
be employed. However, the details will all be omitted. It is known that the 
equations admit one solution for they include as a special case the symmetrical 
orbits. The detailed discussion shows that they admit but one solution. There- 
fore the periodic orbits of the type in question are all symmetrical. 

4. Practical construction op the periodic solutions. 

The demonstrations of the existence of the periodic solutions carry with them 
possible, but very laborious, methods of constructing the series. A direct and 
convenient method of computing the coefficients will now be developed. 

It has been proved that the periodic solutions are expressible in the form 

CO ao 

p> = 2^i^iWi\ w = 2 w.m'; 

(21) *=' '=' 

^' = 2|?;Wi', w' = y^.w'.m\ 

The derivatives are expressible in the same form, and all the series converge for 
all values of m whose moduli are sufficiently small. It has been shown that the 
coefficients of (21) can be uniquely determined so that they shall fulfil the three 
following conditions : 

(a) The series satisfy the differential equation (9) ; and, because these condi- 
tions are fulfilled uniformly over a finite interval in m, they satisfy the differ- 
ential equations identically in m. 

(6) The series are periodic ; and, because the periodicity holds for a continu- 
ous range of values of m, each p^, w., p\, iw^ is separately periodic. 

* Compare Les Meihodis Nouvelles de la 3Iecanique Celeste, vol. 1, p. 88. 
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(c) Since the orbits are symmetrical the origin of time may be chosen so that 
dp I (It = w = dp' jdT = w' = at t = . Since this is true for all values of 
m whose moduli are sufficiently small, each dp./dr, w., dp'./dr and ii)[ sepa- 
rately vanishes at t = . 

These conditions will be referred to simply as (a), (b) and (c). 

Equations (21) are to be substituted in (9) and the results rearranged as 
power series in m . It is most convenient to make these preliminary expansions 
as far as they will be required before the integration is taken up. In the induc- 
tion which follows we shall need the terms in p and w to the fourth power, and 
in p' and w' to the sixth. Up to this order we have from (9) and (21) 

+ (l+;,)(l+«.+ |^-y=[l]4-[2]m+[l + 2 J^-f^,]m^ 

[^(d%i3, dw, dw„ \ /dw„\^ 1 , 

'-ldr+W+P^~di + P^) + [~dr)^-P^-^P4"'+-- 

+ [2^^ - Spl + 4^3 + 2pjm^ + . . ., 

+ "~r (TT/?? [ 1 + 3 cos 2 ( T + «; - «,' )] = "p { [ 1 + 3 cos 2t] 

+ [i^aC 1 + 3 CO® 2t) — 6w^ sin 2t] w^ H } , 

fl + pY 
(22) -I- |- Jf^m^T? -^-- ^-/^ [3 cos (t -)- w — w' ) + 5 cos 3 (t + w — w' )] 

= fJIfj m*?7 { [ 3 cos T -f 5 cos 3t] + • • • } , 

-(i-^)S=-K?]--K?]--[-*55'+4?]-"— • 

-I'^Hr^'^dT+'^'dTdrl'^ +••■' 

(1 + p) 
- ^m^V 7f':^rpy sin 2 ( T + w - w') = - Iwi^t; { [ sin 2t] 



-f [^2 sin 2t + 2w^ cos 2T]m^ + 

( 1 + P )■ . 
— -IJ!/, 7)1*7} pr-— — ^.T4 [sin (t-j-w — -w')-f5sin3(T + M — w')] 

= — f Jf[ rjm* { [sin t + 5 sin 3t] + 



• . • 1. 
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and corresponding equations from the last two of (9) which need not be written 
down. It is a simple matter to write out the expansions as far as is desired. 
In the computations given in this paper the series were carried out so as to 
include all terms in m'". It will be shown that after these series are written out 
there is no need of further reference to the differential equations. 

It has been proved that there are no terms of the first degree in m. This 
result will be accepted though it would have come out in the construction if it 
were not already known. From (22) and condition (a) it follows that the equa- 
tions which define the coefficients of m^ in the solution are 

The general solutions of these equations are 

^2 = Ji? + 2c(^' -f cf cos T + 6f sin T — 97 cos 2t, 

w^ = cf _ (3cf' + j; ) T - 2cf sin t + 2cf cos t + J^i; sin 2t, 



(23) 



^2 = C, T+C2 , 

, /(2) , ,(2) 



"I ' 4^'' ^fi '^f-' '^'x •> '^T^ 1 '^'i^ ^^^ ''f^ ^'■^ *^® constants of integration 
which are so far arbitrary. 

It follows from condition (6) that 



3cf + »? = o, c;® = c;'^ = o 



and from condition (c) that 

c 



f = cf = cf =0. 



Therefore equations (23) become 

^^2 = — i'? + ^2 cos T — 17 cos 2t, ^2 = cf\ 
w^ = — 2c^/> sin T + J5I-1; sin 2t, Wj = , 

where d^'' and c'.^ are so far arbitrary. 

The equations which determine the coefficients of m^ are found from (22) to 
be 

-^^Bp,-2-^'=2-^ + 6p,= -v + 2cf cos r - ^v cos 2t, 



552 F. A. MOULTON : THE PKOBLEM OF THREE BODIES [October 



dp^ dp. 



d^ "' 

dj" - • 

The general solutions of these equations are 

Pg == 2c(,^^ — V + <^T '^^ ''' + ^T ^i" ''' ~ 6 '' ^'^^^ ^"^ ~ ^T"^ ^^° ''"' 

w^ = cf + ( - 3cf + 27;)t - 2cf' sin t + 2cf cos t + i^t^ sin 2t - 2cf t cos t, 

p^ = c^ r + c^ , 

/ '(3) , /(3) 

■«'3 = C3 'r+C^ ' 



where c'-^'' , • • • , c^ are so far arbitrary. 
From conditions (6) and (c) we have 



cf = cf = cf = cf = cf = cf ^ = , 
-3cf' + 27; = 0, 



while d.p and c^ ' so far arbitrary. 
Hence we have 

P2= — 1^ — '?cos2t, 

wjj = Jj-7? sin 2t , 

2).^ = Jt; + 0*2^' cos t — I?? cos 2t, 

TOj = — 2c<2'' sin T + Jg^-T; sin 2t, 

1^2 and zOj are now completely determined. 

It is necessary to go one step further in order to make the induction to the 
general formulas. The equations which determine the terms of the fourth order 
are, from (22), 



P2 


= 


cT\ 


w'^ 


= 


0, 


P'z 


= 


cf\ 


^3 


= 


0; 



d^p 
di^ 



^-3p,-2-^^ = 2^ + 2;,,^-3p| + 6p3 + (^^^^j +8p, 

+ iiiP2+^VP2 COS 2t— 3??W2 sin 2t+^M^7] cos r+^^-lf^rj cos 3t 
= [ - -qW + 1^] + [2cf + |i!f,7?] cos T - |7;[H- 27/] cos 2t 
+ J/-J!fj7? cos 3t + 11?;^ cos 4t, 
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(26) -^ + 2 ^ = - i,, -^ - 2 ^ - 2 -^^- ^-^- - ^vp, sm 2t 

— 3j?m>2 cos 2t — I Jfjij sin t — ^^-M^t) sin 3t 
= [2c(/^— I Jfj^] sin T— f i?(7 +57) sin 2t— i/Jfi?? sin Zt—X^t^ sin 4t, 

Integrating the second equation once and substituting the expression for 
dwjdr in the first, it is found that 

(27) -£i+P.= [2cf - W-'?' + f^] + [-2c^'> + J^J!/,^] COST 

+ |-?7 ( 5 — 3;? ) cos 2t + -Y-J/j J? cos 3t + i/V cos 4t. 

In order that the solution of this equation may be periodic [condition (5)] the 
coefficient of cos t must be put equal to zero, or 

(28) cf = l|Jf.^. 

Thus the undetermined coefficient which remained in p^ and w^ is uniquely 
determined by the condition that p^ shall be periodic. 

After 0^2*' has been determined by (28) the solution of (27) is periodic. When 
the condition has been imposed that dp^jdr shall equal zero at t=0, the 
expression for p^ will carry, in addition to <f^\ an undetermined coefficient <f^^ 
with cos T . When this value of p^ is substituted in the integral of the second 
equation of (26) it takes the form 

^^=[_3ct" + WV-3^] + [3co8T+---. 

In order that w^ may be periodic the constant part of the right member of this 
equation must vanish, which uniquely determines <^*K Then p^ and w^ carry a 
single undetermined constant which is determined in the next step. 
Referring to the last two equations of (26) we find 

X = ¥Tt' + <"t + cf, < = cf * + cf . 

By conditions (b) and (c) we have 
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Therefore p'^, equation (24), is zero. Similarly p'^ can be shown to be zero. 
Consider the equations which define pg and w'^. They are 



2.{E + Mf~ 'i{E'+'Mf 



2 — ^i'* — o / S7 , i/v o / S7 _L u\^ ^''^ ^'^' 



d^w: dEMm 



' sin2T, 



dr' ~2{E+My 

where p'^ = Cj'*. In order that p'^ may be periodic the constant part of the 
right member of the first of these equations must vanish. Therefore 

(29) P, = c - 



2 



6(£: + iIf)'' 
Then it is found, after imposing condition (c), that 

We shall now develop general formulas which will obviate the necessity of 
integrating the equations at each step. Suppose p^, ■ 
P'zi ■ ■ ■' Pi-i' *''2' ■ ■ ■' "^i-i h^-ve been computed. Suppose 



(31) 



p. = a'^"> 4- aj-') cos t + a'^"> cos 2t + • • • + a9'' cosjt, 

w^, = /3</> sin T 4- /S^-*' sin 2t -\ [- ^yHinJT, (i = 2, ■••, 

p'^. = a'^^-'"' + a;*-''^ cos t + a^"^ cos 2t -| 1- a'/-'^ cos ( j — 4 ) t. 



w 



K = /S;^^^ sin T + ;e;"' sin 2t + . • . + ^;i^^ sin (j - 4)t. 



Suppose the coefficients are all uniquely determined except a'j'-^>, /S'^*"'', which 
we shall suppose carry an undetermined constant linearly, and a'^ and a'^ , 
which we shall suppose are so far arbitrary. The problem is to find the coef- 
ficients in the expressions for p^, w^, p'. and w[, and to determine a'/'''^ ^i~^\ 
a^ and a^ 

Consider the differential equations for the determination of p^, id., p\ and w\ . 
The left members of the differential equations are of the same form for all 
values of i beyond 1 , for example, as given in (26). The right members are 
polynomials in /»,,••• , p._^ , w, , ■ ■ • , w._^ ,p[,---, p'._^ ,w[,---,w'._^, and linear 
in cosjT and sin jV. In any term which is independent of cosjr and sin ;t the sum 
of the products of the exponents and subscripts cannot exceed i. In any term 
which involves cosjV or sin jV the sum of the products of the exponents and 
subscripts cannot exceed i — j. It follows from the way in which j), w,p', w' 
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and T enter in the differential equations (9) that the right members of the first 
and third equations which determine Pi,^.,]}'. and w'. are even in the dp./dr, 
w., dp'.jdr, w'. and t taken together, and that the second and fourth equations are 
odd in the same quantities. Hence it follows from the forms of equations (31) that 
when the values are substituted for the p., w., p. and w'. , the right members of 
the first and third equations will be sums of cosines of multiples of t, and the 
right members of the second and fourth equations sums of sines of multiples of t . 
It is necessary to find how P;_i, W;_, , and p!_, enter in the right members of 
the differential equations which determine p^, ?«., p'.^ and w\. It is seen in the 
special case (26), and it can be shown in general from equations (9), that the 
first equation contains the terms 6/?;_j and ^dio._^jdr, the second — 1dp._^jdT, 
and the third Zp\_^. These are the only terms that contain coefficients which 
are not completely determined. Hence the equations which determine p^ , w. , p\ , 
and w'. have the form 

t 

-^ - 3p, - 2 -^ = Af + ( 6<-» + 2;8(;-') + A-\) cos r 

+ J.'/' cos 2t + 1- A^f cos iV, 

(32) ^^ + 2 -jj^^ = (2af-'' + Bf) sii> T + 5« sin 2t + 1- Bf sin ir , 

5r=(3<'- +^;"->) + ^f COST+ ...+^f:',cos(^-4)T, 

^^Ji = B[^'^ sin T + . . . + 5:^*, sin (* - 4 ) T, 

where the A'-p, B^, A'}' , and B'. '' are known numbers depending upon the 

coefficients of ^^^ ' • ' W,-i' '"'2' ' ' '■ '"'.-i W2' • ' '» Pi-i' ^2' •••^Wj-i- The 
constants «i'~'\ yS'^'""'' and a,', '~^' are as yet undetermined. 

The solutions of (32) are to be periodic, and dpjdr, w^, dp'Jdr, and w! must 
vanish at t = . Therefore they must have the form 



(33) 



p^ = c^p -|- a'j*' cos T + <^p cos 2t + • • • + a^ cos tr, 

w. = ^p sin T + ^f sin 2t + • ■ • + /3'^'' sin it, 

p)^ = a'p + a.'^^ cos T + a'p^ cos 2t + • • • + a'^2^ cos ( i — 4 )t, 

w\ = /3;''' sin T + /Sf ^ sin 2t + . • ■ + ^[^^^ sin (^ - 4)t. 



The dp. /^p, OLp , and /SJ '' are to be determined in terms of the 
A'-!-', Bp, Ap', B'. ' , and the new constants of integration. 

Consider first the third and fourth equations of (32). In order that p'. may 
be periodic we must impose the condition 

Trans. Am. Math. Soc. 37 
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Then imposing conditions ( 6 ) and ( c ) we have 

a'J'^ = constant undetermined until computation of p^_^^^ , 



(34) 






('■) 



0- = l,.--,i-4), 



Now consider the first two equations of (32). Integrating the second we 
have 

(35) ^ = cf - 2p. - (2a<;'"'^ + Bf) cos t - ^B^^ cos 2t v Rl" cos it . 

Substituting this expression for dwjdr in the first of (32), we find 



^Pi 



--jii+p^ = [A<-p + 2c(/>] + [Af> - 2Rp + 2a(;-'^ + 2/3(*-i'] cos r 



cZt^ 



(36) 



+ 
+ 



^2 2 2 






cos 2t + • • • + 



[^''-'^'] 



cos^'t + • • 



cos IT. 



In order that the solution of this equation may be periodic [condition (6)] we 
must impose the condition 



(37) 



A'{> - 2B^p + 2a</-i' + 2/S^;-') = . 



Then the most general solution of (36) satisfying condition (c) Idp./dT = at 
T=0] is 

p. = [ ^(;' + 2cW] + dp cos T - ^-2^:^! ) ''"^ ^ 

(38) 

— *- .i . ; ?-T^-^ cos ;t — ^ ~- •7-.2 — i^N cos it . 

j{f-l) -^ ^(^2-l) 



Substituting this value of p^ in (35), we get 

^^' =- [2^(,'->+3cV)]- [^V''+2<+2a(/-'>] cos t+ '—2-/22—1 )^^o« 2t+. • • 
[2;^y)-(/+3)i?«] [2i^f- (^^ + 3)^'/'] 

+ Jlf^TJ COS^T+...+ i(j2_l) cos IT. 
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In order that the solution of this equation may be periodic we must impose the 
condition 

(39) 2^« + 3c« = 0, 

which determines c*-^^ uniquely. Substituting in the first term in the right 
member of (38), we find a^'^ = J.^*^ + ^c^/^ = — iAl'\ Then the most gen- 
eral solution of the equation for w. , vanishing at t = , is 

w, = - [^;) + 2a(;) + 2af-')] sin t + !- ^-j^ '—-'-sin 2t + . • • 

(40) 

+ /(/-I) "'^•^^+ • • • + ^~~^ Ll) "" ''■ 

In the notation adopted in (83) the coefficient of sin t in equation (40) is ^p. 
Since the form is general we have 

(41) /3(;-^) = - [ ^/-i) + 24'-') + 2oi^r''' ] • 

This equation used with (37) determines a\*-'^'> and ^^-'^'> in terms of known 
quantities. Hence the general expressions for the coefficients of (33) are 

<'=-i^i'' (i = 2, ..., ^), 

«(*-'> = 1 [Ap - 2i5W] - [5('-') + 2a(*-2)] , 
/3(^-i)= - [Ap - 25W] ^. [^(i-i) 4. 2a(/-2)] , 

ri^f-2».')3 



iO'^-i) 



(42) 



^■■' = 


[2j4"-(/ 


+ 3 


)^^ 


/(/- 


-1) 




,,fi-2) _ 
'0 


- iA':- 






af = 


i^,« 

z^^- 






^f = 









(i=2, ••■,»), 

(i = 6, ■••,00), 

(i = l,---,i-4), 



The simplicity of the computation from these formulas and the expansions 
(22) will be illustrated by the computation of a^f^ and ^'•9. 
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Table I. 



Term from (22). 


Constant. 


COST 


cos2r 


cos3r 


cos4r 


+ ^1? 





- 2^^' 


+ ¥'? 








+ =^.S- 


--W 





-\w 





- ^'^v' 


+ 6p3 


+ 2,, 

+ Wv' 


+ 6af > 















+ 3p, 
-3pl 
+ i'?A(l + 3cos2T) 
— SrjWj sin 2t 
+ f »;Jfi(3 cos T+ 5 cos 3t) 


-h 

-w 

-iW 











-3, 















^j« 


+ h-Wv' 


+ h^. 


-%v-h' 


+ -'^vM, 



Table II. 



Term from (22). 


sin r 


sin2r 


sin 3r 


sin4T 


P^ dJ 





-kW 





--VV 


odp, 
-^dr 


+ 2a<^' 


--VH 








dr dr 











--VV 


— hP2 sill 2t 





+ lrf 





+ 1';' 


— SrfW^ COS 2t 











-fIV 


— f »;i/i(sin T + 5 sin 3t) 


-hM, 





--V-^M, 




1 5 3„2 


Bf 


-%r,M, 


--VS-I'?' 


-J/^if, 



From these values of A'^f' and ^.*' ( ^?' and a^^^ are zero) it is found by equa- 
tions (42) that 



2 



^,(4) '- _L 3 3 1 „2 



(43) 



f= + lf^Jf„ 



7(-l)_ 



1 0„ I 2-1! 



As many succeeding dP and /3'.'' as may be desired may be computed in precisely 



,(4) _ _ 2 6 
"s — "6 4 



V3I„ 



"4 — 8'/ ' 
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the same way. The operations are reduced to systematic routine, and the work 
consists almost wholly in multiplying fractions in constructing the tables, and in 
adding fractions in getting the A'-P, MP, J.'.'*' and 5'.'*'. 

The differential equations admit two integrals of which no use has been made. 
They may be used to derive new formulas for the computation of the dp, ^P, 
a'. ' , and fi[!' , or to verify the results obtained by (42). The principal applica- 
tions occur when the mass of Ifis infinitesimal, or at the most so small that its 
effect on the motion of S around the center of gravity of jE and M is insensible. 
Consequently, to abbreviate the discussion we shall assume that S moves in an 
undisturbed circle around the center of gravity of JiJ and M . Then the integral 
of areas disappears and the vis viva integral becomes Jacobi's integral: 



(44) - 2»^' 



H«7 I ^jl-2^,{l+p)cos{r + w) + (^^yi+py 



a 
(1 + p) cos ( T + w ) 



= C, 



where C is the Jacobian constant of integration. 

In the periodic solution p and w are developed as power series in m, and 
hence we may write the integral for this solution in the form 

„ ^ ( dp, dw. \ 



= C. 



Since this integral converges for all values of m sufficiently small, it follows that 
each F^ is separately constant. In F^ the sum of the products of the exponents 
and subscripts of the factors of any term independent of cos_;t cannot exceed i ; 
and in any term involving cos jr the sum cannot exceed i —j. The function JF 
is even in dp./dr, w. and t and therefore each F^ has the form 

(45) F. = <y\l'> + 7<j'' cos T + • • • + y'-P cos jt + • • • + 7'/' cos It. 
Since every F^ is identically constant in t it follows that 

(46) 7f=0 (i = i,..., ,x,;i=l,.--,0. 

The fulfillment of these relations, which are functions of dp, ■ • •, a!f, jS^p, ■ ■ -, 0f 
(J — , ■ • ■ , i ) , serves as a complete check on the expansions and the computa- 
tions. 
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The integral is applied by first expanding the separate terms of (44) as power 
series in m. Thus we have the expansions 



m 



.^ ( 1 +pf= —m^ — ip^m*-^- 



(47) _ ^A1±^J = - 2 - 4m - 2 [ 1 - ^,] «j^ - 2 [ 2^, + ^3 ] m' 

— <i[p^—pl + 2^3 + ^J m* + . • • , 

^(1 +/>)cos (t + w) I = — i77[l + 3 cos 2t] m^ 

+ [— jyJjC -'■ + *^ cos 2T) + 3??t02sin 2t— ijf,7?(3 cos t+5 cos 3T)]m*H . 

It can be shown from these functions that in general F^ contains the terms 
4cp^, 2dw./dT, 4pj_^, and terms with lower subscripts. Suppose the coefficients 
of p^, ■ ■ ■, pj_2i '^\i " '■> '^'.■-2 ^^^ entirely known, and that the coefficients of 
Pf_^, and w^_j are all known except a^*"'' and /S['~". Let in general the known 
part of 7^' be C-P. Then the conditions (46) become 

4aW + 2/3';' + 4af-i> + C« = ( ,; ^ 2, • • • , oo ) . 

^ ^^ ia!p + 2j8f+Cf = = 2, ■..,{), 

From equations (42) it is found that 

4ctf + 2ySf + 4a(;-'' + 2Bf = (i=^2, •••,00), 

4af + 27/SW + ?^)=0 (i = 2, ■■:,»)• 

The comparison of these equations with (48) shows that 

2 

(49) C«-^.^<«=0 (i = 2,---, oo;j = l, ■.■,^)■ 

When the B^''} have been eliminated from equations (42) by means of these rela- 
tions, the oc'?' and /3'*' are expressed in terms of the A'9 and C!*' by the equations 
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,«) _ 






(60) 



f-i- 



(j = 2, 



(i = 2, •■•,0, 



/3<« = 



2j(/-l) 



Therefore in case the orbit of S is a circle around the center of gravity of £J 
and If the second differential equation, which involves the second derivative 
of w, may be entirely neglected in the determination of the coefficients of the 
periodic solution, and Jacobi's integral used in its place.* The more conven- 
ient process is to be employed. In every case there are more terms in the table 
for the C-P than for the £'•9, but a large proportion of the former are the same 
as those occurring in the table for the AJ^ except for numerical factors. For 
this reason the a!'!-' and ^^9 are more conveniently found from equations (50) 
than they are from equations (42), but then the advantage of having the inte- 
gral as a check is lost. In the computations given in this paper the coefficients 
were obtained by (42), and the integral was used as a check. 

The use of the integral will be illustrated in computing the coefficients of m* . 

Table III. 



Term from (47). 


Constant. 


COS " 


cos2r 


COST 


cos4- 


Ht-J 


+ 27,' 











1 


+ip, 


+ 4a<*' 


+ 4:a[» 


+ -H" 


+ M'> 


+ 4a[*' i 


-pi 


1 9„2 

— 3 6'? 





-w 





-w 1 


■^"-t- 





+ 2^;*' 


+ 4^8^*' 


+ OyS?' 


+ 8^8"' 


+ m 


+ Wrf 











+ Wv\ 




--VV 





1 1„2 





-i^rf 1 


-VP2 


+ H 





+ v' 








+ *Ps 


+ h 


4af> 


--h*-v 





1 


— 3r]p^ cos 2t 


+ w 





+ iv' 





+ N \ 


+ 3r)W^ sin 2r 


-4- ^S^n' 











-iW 


— iil/,»? (3 COS T 4- 5 cos 3t) 















cr 


—kS-W 


1 5 3„2 



* Hill used the integral to eliminate the most troublesome expression arising in the differ- 
ential equations in his method of treatment. Loc. cit., p. 132. 
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With these values of the C-f^ and the Af^ of table I, equations (50) lead to 
the results of (43) ; or, when the values from (43) are substituted in table III, 
the relations (46) are satisfied. 

Suppose the a^p and /3j*^ are computed by (42) and that the integral is used 
as a check. It follows from (48) and (49) that at the ith step the check does 
not cover the A'^^ . That is, an error might exist in any A'^P without its being 
revealed by the check formula until the ( i + 1 )st step. However, there is the 
check for the A'-p and C','-' that the sums of the coefficients in every line of 
their tables must equal the values of the functions in the first columns at 
■T = 0, and the sums of the A'-P and C-P (with respect to j) must equal the 
sums of their respective first columns at t = . 

6. Applications to the lunar theory. 

The differential equations and their solutions have been developed in such a 
form as to be immediately applicable to the Lunar Theory, provided we assume 
that its motion is periodic of the type under consideration. It is only necessary 
to find the order of magnitude of the parameters involved so as to know how 
far the developments must be carried to meet the needs of observers. Observa- 
tions have given for the Lunar Theory the following approximate numerical 
relations : 

n 1 

n = 13.4w. , whence m = , = ^-^ ,- , 

71 — n 12.4 

B= SOM, 

a i 

(51) S = 330,000^, whence v = "o b^ rr = 

O -)- jQ + -M , 



1 + 



300,000 



- — ^1-- whence M — ]V — 9- M — JV = l M — JSf = -^^ 

Consider first the expressions for -p and w . The term of lowest degree in 
J) is, by (29), p = p\ m^ = EMN,m*l<o ( E + Mf . It follows from (51) that 
p'^m"^ = 1/80,000,000. That is, since r' = a'(l +/)'), the average distance of 
the sun from the center of gravity of the earth and moon is about 1.2 miles 
greater than it would be if it were not disturbed by the moon. The first 
variable terms in p and v) are, by (30), 

p, m" = + 1 j-^-^j^Y ''^ ' "*« m" = - f jE-:^Mj *'" ■ 

Therefore, with respect to the center of gravity of the earth -and moon as an 

origin, the sun describes in a period of one half a synodical month in the retro- 

grade direction a nearly circular orbit of radius 95 feet around a point whose 

distance is approximately a {1 -\- p[m* -\- p'.^m' + p'^^m^^ , and which revolves 

with the mean angular motion of the sun. 
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The coordinates of the sun, p and w' , enter in the right members of the first 
two equations of (9). It has been shown how slightly they vary because of 
perturbations by the moon. The reciprocal action of these perturbations upon 
the moon are still smaller, because where they appear in the differential equations 
they are multiplied by m^ . If we regard MN^I^E + M) as of the second order 
in m , the perturbations of the moon by the sun due to the fact that the moon has 
disturbed the sun from a strictly circular orbit are of the 8th order in m, and 
their numerical coefficients are so small that they are absolutely inappreciable 
observationally. That is, no sensible error is committed in computing the 
moon's coordinates by assuming that the sun goes around the center of gravity 
of the earth and moon in an undisturbed orbit. 

It is seen from (51) that the constant rj in the Lunar Theory is very nearly 
equal to unity. Since the higher terms are much simplified by replacing it by 
unity while no sensible error is committed, we shall retain it explicitly only up to 
terms of order 4 in m. By applying equations (42) where the ^^*^and B'P are 
defined in (32), it has been found that* 

^ = - 1^ g + cos 2tJ vm' + 1^ 3 + jg Jf, cos T - g cos 2tJ vm' + ^^^ r, - ^ 

3 2 25 

+ oo (45?; — 1 )ilfj cos T 4- Q ( 3»; — 5 ) cos 2t — g^ Jfj cos 3t 

3 , 1 4 r 181 3145 ,^ 19 „ 95 ,^ ^ 

— g r? cos 4t 7;m^ + jg- ^ -)- ^^^ M^ cos t — ^-^^ ^^^ 2t -)- ^^3 ^i cos 3t 

, 1 , r/ 3175 225 ,^, 3 , \ 56831 ,^ 
^"J "" + LV 32"8l + 512 ^? - 16 ^^) + I02¥ ^1 ^«^ ^ 



32 -3 -5 



_4897 
16*^9^25 



/ 10933 225 ,^, 6 \ „ 5053 ,^ „ / 

- ( 256T8T + 512 ^> + 8 ^^j ^«^ ^^ + 256^3 ^^ "'^ ^^ " ( 

7 ,^\ , 155 ,^ ^ 59 ^ "1 „ r/ 245 105 ,^, 
+ 32 ^,) cos 4t - g^-2 ^1 cos 5t - g^g cos 6tJ m^ + |^(^ g^^Y + 32 ^> 

^ .A / 910^469,^ 106 ,^\ / 911,671 8985,^, 

+ 8 ^^ j + ( 1024; 3 ^^ + 128 ^^3 j «o« - + ( 5T2T243T5 " 1024 ^' 

5 , \ „ 285,283 ,^ „ / 130,723 525 ,^, 

- g M,j cos 2t + ^^gg^ 3 Jf, cos 3t - (^ 32^.^rri25 " r024 ^^ 

49 ,A . 625 ,^ , 51853 , 1 , 

+ 32-3^5 ^^j '^"^ *^ + r024^3 ^^ "^^ ^^ " 512-357 °°' ^^J "" 



*Mr. W. D. MaoMillan has duplicated the entire computation as a check. Moreover, the 
integral and other check formulas have been exactly satisfied, and it is not believed that there 
can be any errors in these series. 



564 F. R. MOULTON : THE PEOBLEAI OF THREE BODIES [October 

r/ 2,228,365 165^5 53 \ /228. 305,129 1125 

■^ |_ V8192^43 + ^192^' ^' + 31 V + V ~T6384 • 9 ' ' + "25T • « 

225 ,^\ / 8,581,859 203,593 ,^, 1681 ,^\ 

+ 64 ^^sj cos T + (^ gj-^29T-2T - ^09r ^^ + -256" ^^j '^"^ 2^ 

1,145,073 ,^ 245 \ „ / 8,398,553 4245 



/ 1,145,073 ,^ 245 \ „ / 
+ ( W96V3-^' - 10-24^3) cos 3t - ( 



If 2 



+ 64-9 



.256 -81 -625 2048 ■ 
324] 



1 ..\ . / 943,835 ,^ 147 ,^\ 

725 "^^ j *^°' ^^ + V 16384 -9 ^' " r0^4 ^^^ j '^^^ ^^ 

/ 698,591 425 ,^, 49 , \ ^ 31885 ,^ 

- (3-2^9725"^ + 8192^' + 256 ^^j "''' ^^ - 16384 -^^^ *^"^ ^^ 



„ -1 , r/ 13,799 5,765,369^^, 5 ^ \ 



/ 4,337,803,547 ,, 7575 ,^ ,^ 64,025 ,, 166,875 ,^ \ 
+ i "-6553-6 -9- ^' + ^5T ^^' ^^^ + 2048^3 ^a + "32768" ^^' j ''"^ " 

/ 1,924,392J99 _ 96:8^3 346,881 \ 

^ V512 -2187 125 512-9 ' ^512 -9 -5 V 

^1,018,308, 773 14,445 3045 10,125 \ 

+ 1^ 16384 •27-5 ' ~ 8192 ' ^2 + 4096 ^ + 32768 V *'°'' 

/ 1,092,030,283 _ ^348^655 _?^^L m~\ 

+ V r024727"-3r25 -^7 "~ 16384^9^5 ' ~ 64727^25 ^ "^^ 

/ 264,181,703 ,, 3045,^,^ 217 ,^ \ 
+ [ 65536-T2TT6 ^^' + 8l92 ^' ^^^ " 4096 '^^^ j cos 5r 

/ 1,597,008,893_ 62,915 3271 \ 

+ V ~ 5T2^27"1.26~343 "^ 4096" 3 -7 ' ~ SlTlTS V ^°^ 

5057 ,^ ^ 5,778,113 „ 1 , 

- 3'2768 ^' '^^^ ^^ - 8192T3"5T49 ""' 8^ J "^^ + ' " • ' 

V = ( 1 + OT ) T + -o^ sin 2t ?;ot.2 + 7^ J/j sin t + -^ sin 2t »;m'' 

[3 1 15 201 1 

- jg(45r; + 7)il/j sin t— jg(9»;-41) sin 2T+g2^¥j sin 3t + ^ .g sin 4t hm* 

r 4151 ,^ . 89 . ^ 25 ,, . „ 2177 . '1 , 

"^ ~ "128 ' ^'"^ '^ "^ 4^7 ^"^ ~ 16 ' *"' "*" 64" 3 ■ 5 *^"^ 
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r 71,449,^ . / 62,911 1125,^, 25 _\ 

+ [ ~ ^512 ^^ '''' " - i2648T81 " 1024 "^x "32^0 ''° ^" 

413 ,^ . ^ / 23,737 63 ,^\ . ^ 75 ,, 

- 32 M, sin 3t + (^ 32^^-25 + 256 ^^ ) '"' *^ + 128 "^^ ^'" ^^ 

^ ^^1^ • «1 6 , r /2'^90'245,^ , 105 \ . /- 578,077 

+ 2048 .3 sm 6tJ m^ + [^- {-j^^^;^ M, + -gj Jf,j sm r + (128T243-5 

34,155 5 \ . 100^995 . „ , / 1,339,541 

+ ^048 ^^ + 4^0 '''' ^" ~ -2048"^^^ ^^'^ ^^ + ( Iir-lTTiS 

1725 ^„ , 77 ,,\ . , 919 ^^ . ^ 4289 1 

- 2048''^^ + 32 :3 .5"^^ j '''' ^^ - rm^^ '''' ^^ + i6^3~5 -1 '''' ^'^J'^' 

r / 90,908,791 1125 ;^. ;^ _^ 135 \ . / 54,988,283 

+ L ~ V ^192 • 3 " ^^ + T28 "^' ^^ + T6 ^3 j s'n ^ - ( 



256 -729 • 25 



811,967 _ 45,721 ^^\ . „ /7,378,645 ^^ 385 \ 
8r9T-^^^ + 204Fr3 ^0 ''" ~ V~409trT^' " 512T3 ''^V '^° ^' 

/ 652,081,58 3 10,12 5 _91,069 ^\ . 1 / 133,763 ,, 

"^ V 8192 • 81 • 625 ~ 4096 ^ ^ 512 • 9 ■ 25 ^^ j ^m 4t - ^ " 8192 ~ ^^ 

399 ,^\ . ^ / 338,729 1675 ,^„ 693 ,^\ 

- 512:5 ^3 j «>° 5t + (^ 256-5T49 + r6384 ^J + 2048 ^^0 '"' ^' 

279,055 ^ . „ , 207,155 . „ 1 « f/ 15,772,500,205 
+ 8192.49^' ^^° ^" + m0T-2T3«'"8Tjm« + ^(^-^32^^-g-27--^^x 

^'^^^^^MM ^'^^^'^^ ^r l'^3,625 \ . / 1,601,782.811 

- -256 ^^^^^ - 1-0X4:3 ^^3 - T6384 ^' j ^^° ^ + (" 256 .2187^125 

, 4198 182,881 „ \ . „ / 537,179,239 ,^ 6075 

+ -^9-^^-5l2T3--5^0"° '""^ (- 4()96T27:5^^x + 4096^^^^^ 

_ - 9205 14,625 \ / 3,297,305,351 

2048 • 3 3 16384 'J ^"' + V ~ 512 • 81^3125: 7 

5,397,281 J^7,537 ^\., / 22,567,883^ 

■^ "32768 ■ 9" ' + 128 • 27 125 ^^' J ^'° ^^ + \^ - ~S27W-W~ ^^1 

22^^ uw^'^^'^ Af\ ■ ^ ^( 1^673,061,443 195,865 ,^, 

- 4096 ^^ + 2048^;^'"^" + (,64^81:7125:343 " 16384:7 ^' 

5659 ^^\ . ^ 6857 ,^ . „ 

+ 512: 375 '^^) '"' ^^ - 81927 37-7 ^' '"' ^" 

70,605,533 . „ T 
+ 32768:3:5.49''"^^]"^ +•••• 
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If in these series r] is put equal to unity the terms which are independent of 
3f^ , Jfj > ^^^ ^3 correspond to those found by Hill in his researches on the 
Lunar Theory.* It is not easy to compare this work with Hill's for he used rect- 
angular coordinates, and in the final series a different parameter. Moreover, 
his linear constant was different from the a defined in (4) and used here. With 
7j=l the terms involving Jfj, M^, and M^ correspond to those found by BEOWN.f 
The difficulty of comparing with his results is the same as that of comparing with 
Hill's work. 

Delaunat used the reciprocal of the radius vector, the longitude, and the 
latitude as coordinates in his work on the Lunar Theory. Instead of the para- 
meter m he used in its place the ratio of the mean motion of the sun to that of 
the moon. If we call it /i, it is 

n' fi 

a= — . whence m = ;; . 

n 1 — fi 

Delaunat and other investigators in the Lunar Theory have used Fourier 
series whose coefficients are power series in the parameters employed. If we 
rearrange the series above as Fourier series and transform from m to /i , we may 
easily compare with Deladnat's expression for the longitude. In order to 
make the comparison we must, in the series above, put 

The sin jt above is sin jD in DelaUNAy's notation. Perfect agreement is 
found in the coefficients of all terms except the following : 



*American Journal of Mathematics, vol. 1, p. 143. 

t American Journal of Mathematics, vol. 14, pp. 140-160. 
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Coefficient of 



a 

a 
a 



/jt sin T 



fj/ sin T 



(j)M«sin 
fi' sin 2t 
fj,^ sin 2t 






sin 2t 



— ; M sm 6t 
a 

— r«' sin or 
a 

(^)Vsin3r 

( — ; I w* sm dT 

fi,^ sin 4t 

/i® sin 4t 

(j)/.^sin4T 

(^)M=sin4T 

— ^ M° Sin 5t 
a 

—7 «" Sin 5t 
a 

-7 M sm 5t 



Delaunay's value. 



63,106,813 
8192 
10,835,537,15 9 
" "66536 -3 

not computed 

7,596,606,727 
"•" 4096 • 729 25 

8,051,418,161 
■2048-2187026 

not computed 



(( (( 



(I (t 



tt (( 



+ 



123,030,377,303 
8192 -81 -626 

not computed 



3911 

+ 1024 

not computed 



Present computation. 



74,009,397 
~ 8192 
47,823,726,517 

32768-27 
2415 



+ 



128 

6,981,835 

~ 16384"^ 

3,748,175,501 

2048-729-25 

475,28 1,673 

128 • 218r-T25 

3,713,357 

+ 



+ 



8192 3 
86,194,253 
2048-9-5 
21,715,556 
~ 4096 - 9 
1,764 ,511,03 9 
4096-27-5 

385 

+ 



+ 



+ 



612-3 

16384-3 
82,679 ,012,303 
"Sl^^^'Sl^T 
876,528,8 10,3 03 
1024-8r~3125^ 
_ 2,754,403 
" 4096^"9~25 

1,217,563,507 
■32768^7^125 

1481 

-f 



1024 
122,623 
" "Ti9 2^ 
51,974,771 
32768^9" 
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Coefficient of 



\ , ] fj'' sin OT 
\a J 

I —, 1 w' sin OT 
\a J 

/u.' sin 6t 

fjiy' sin 6t 

{ - 7 I u* sill 6t 



Delaunay's value. 



sm OT 



, fjfi sin 7t 



7 /U.' Sin 7 T 



/A* sin 8t 



/u.' sin 8t 



not computed 



4( fci 



&( i.i. 



fc( i; 



11 u 



il. v. 



(1 u 



u u 



u i( 



bb bfc 



Present computation. 



+ 



399 

1309 
+ 4096 
54,143,779 



2048 -3 -5 -49 
138,760,83 8,897 
256^81^25^^43 
7219 
+ 16184 
1,361,621_ 
~ i6384"- 3 ~5^ 
279,055 



+ 



+ 



8192-49 

4,97Jt,9^1 

'8l92~1^9 

207,155 



+ 



+ 131072-3 
34, 667,316,70 3 
32768 - 3 - 5 - 4^ 



The very general agreement shows that Delaunay's terms involving the 
argument ^'D belong strictly to a periodic orbit, and that he carried out the 
computations with remarkable accuracy * considering the complicated operations 
which his method demanded. 

When the expressions for the coordinates are arranged as Fourier series, it 
is found that the coefficients expressed in terms of m converge much more 
rapidly, at least to the ninth power, than when expressed in terms of /a . Hill 
showed t that if a parameter ni is defined by 



m = 



1- 



the convergence is still more improved. This transformation removes the pole 
arising in the expansion of the principal term from | m | = i/6 to | «i' | = 3 . 
In /A the pole is at the distance | /-t | = l/^ i • ^^ ^° '** ' ^® introduced by 

*Beown's results disagree with Delaunay's in the coeiBcients of {a I a' )ii^ sin jD; he did 
not compute higher powers. Loo. cit., p. 155. 
t Loc. cit., p. 141. 
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m = 



6 — 4m + m^ 



the pole of the expression is removed to infinity and the convergence of the 
series for the coefficients of cos 2t and sin 2t, at least out to the 9th power, is 
slightly improved. But these coefficients are given by any of the series with suf- 
ficient accuracy. The slowest convergence is in the coefficient of sin t. In 
Delaunay's parameter this coefficient is 



[- 



15 93 , 6887 , 137,197 , 4,028,333 , 74,009,397 

; ^ 8192 

47,823,726,517 



^ ~ « '^ ~ 128 '^ ~ 512 ~ '^ ~ "1024-3 ''' ^8192 ^ ^ 



32768-27 
In the parameter of this paper the same coefficient is 



7 1« 
J a 



(52) [- 



15 39 , 4151 , 71,449 ^ 2,190,245 , 90,908,791 



m— -VA-?n — cT-ni ew-t^ttt; — W' 



4 128 512 1024-3 8192-3 



15,772,500,205 , 



32768-27 



]a 
a 



When expressed in seconds the coefficient of the term containing the 6 th and 
7th powers of the parameter is — 0".7122 and — 0".3811 in Delaunay's 
series, and — 0".5453 and — 0".2124 in the series given here. While the two 
series represent the same function, the first seven terms in m seem to give it 
more accurately than the first seven terms in fi. 

It is noticed that the numerical coefficients so far as given are all of the same 
sign, and that the ratio of each coefficient to the preceding one seems to tend 
approximately toward 5. If this were the precise ratio, and if it were indefi- 
nitely maintained, the series would be simply the expansion of (1 — 5m)"'. If 
the first few terms succeeding those which have been computed proceed approx- 
imately according to the same law, then the function will be much more accu- 
rately given by them if they are multiplied and divided by 1 — 5m . Thus we 
find by these operations in place of (52) 



(53) 



15 3 , 2089 , 11.571 , 46,775 , 3,298,991 „ 
ym-gm--f -^-2^mH -5Y2- "^ " 1021^ ™ " "Z^M^^ '"^ 

591,082^75 , 

+ ^32768~27" '"' ' 



In applying Hill's method to the parallactic inequalities Brown found* that 



* Loc. cit., p. 149. 
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the divisor 1 — 4??t — ^^m?' • ■ ■ naturally enters. The terms in the series with 
this divisor decrease somewhat less rapidly than those in (53). It is to be noted 
that the radius of convergence is never altered by using any finite polynomial 
as a multiplier and divisor in this way. The process is in reality a sort of 
method of exterpolating which is useful if the coefficients proceed for a time 
according to the assumed law. 

The average ratio of each coefficient to the preceding one for the six terms 
given in (52) is approximately ^^-. Hence we may conveniently write the series 
as — -y-TO[l —^^-m^~^ plus a series which shall make, when added to the ex- 
pansion of this function, the six terms given in (52). We find in place of (52) 

r 1- r. o, n , „ , 4637 , 549,839 , 

(64) 

48,596,2 07 ^ 2,499^42,713 ^ 174, 308,031,761 ^ la 

+ 1024 • 3 •'125 ™ "^ '8192^3"^62F"^ ~ MWS ■ 27 • 3125 "' '^ "'\a' 

To illustrate the effects of these processes we shall reduce these series to 
numbers. It has been found from the observations that 

m = 0.08084893, -^ = 0.00255878. 
a 

With these numbers we get for (52), (53)^ and (54) respectively 

- 80".008 - 33".637 - 9".045 - 3".147 - 1".300 - 0".545 

-0".212... = -127".894, 

- 80".008 - 1".294 + 4".552 + 0".516 - 0".028 - 0".020 + 0".008 • • • 

0:59576535 

= -128".029, 

- 127".382 - 3".882 + 2".021 + 0".969 + 0".231 + 0".024 

-0".001... = -128".020. 

The corresponding values in Delaunay's series including a term of the 7th 
degree in fj, are 

- 74".024 - 34".330 - 11".886 - 4".428 - 1".862 

- 0".712 - 0".381 . . . = - 127".622. 

Brown found * for this coefficient, using the numerical value of m from the 
outset, — 128".069. W^hile this number is slightly larger than one would infer 



'Monthly Notices of R. A. >S., vol. 64 (1904), 534. 
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from (52), (53), or (54), it is -very probably essentially correct. But in the 
computation we put rj = 1, M^m^ = a J a', whereas the correct values are 

. = --1-^ = 0.999997, 

1 + 



M,m^ = 



S' 
E — M a 80.487 a 



' ~ B -{- M a' ~ 82AS1 a" 

With these values of the constants the coefficient (53) is — 124".926. Hill 
has computed f the coefficients of the various arguments of the series with a 
very high degree of accuracy by substituting their approximate values as deter- 
mined by himself and Brown in the differential equations, and giving them such 
differential corrections that they more exactly satisfy the equations. He found 
for the coefficient of sin t, when we change from his value of a fa' to that used 
here, — .124".965. 

6. Applications to Dakwin's periodic orbits. 

In his memoir in Acta Mathematica, vol. 21, pp. 99-242, Darwin con- 
sidered two finite bodies of masses 10 and 1 revolving around their center of 
gravity in circles, and an infinitesimal body moving in the plane of motion of the 
finite bodies and subject only to their attraction. He discovered a large number 
of periodic orbits by interpolation from series of special orbits computed by 
mechanical quadratures. It will be shown that his " Satellites A " and " Planets 
A " belong to the type of orbits treated in this paper. 

To make the series computed in this paper applicable to the satellite orbits 
treated by Darwin it is necessary to put M=(i, E=\, S =1Q, a' = 1. 
Then, when the infinitesimal body moves around E, we have from (4) and (6), 

S 10 

''~:r+^~ii' 

(55) 

Darwin defined his orbits by the value of the Jacobian constant, and their 
periods were found from the detailed computations. The Jacobian integral may 
be written 

F= F^ -f F^m -f F^m'' + . - . = C. 

For an arbitraiy value of m sufficiently small, as m^, the value of C can be com- 
puted. Let it be C„; then 

* Astronomical Journal, vol. 15 (189 ), p. 137. 
Trans. Am. Math. Soc. 38 
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Let the value of C of the required orbit be C = C^ + A C, and the correspond- 
ing value of m be m = wij + Am . Substituting these values in the integral and 
inverting the series, we have Am expressed in terms of A C, and therefore m 
expressed in terms of C. With this value of m the coordinates can be com- 
puted from the series which have been given above. 

While there is no difficulty in carrying out the details of the steps just de- 
scribed, it will be sufficient for the comparison to take the period as found by 
the computation. We shall compare with "Satellite A" (Jacobian constant 
= 40.6), loc. cit. p. 199. The synodic period was found to be 61° 23' = 61° .383, 
where the period of the finite bodies is 360°. Therefore 

61-383 1 ^-.^^c, 

"^ = "^0^= 5:8648= ^-^^^^^^ 



a 

a 



- = (n)*(ri-»)'— 



The m is more than twice as large as that occurring in the Lunar Theory, while 
the a/a' is more than 40 times as great as a/a' for the moon and sun. 

Putting M.m^-' ■= {^a/a )' in the series, and using these values of the constants 
we find 

r = 0.12427 -f 0.00652 cos t - 0.00420 cos 2t -f 0.00004 cos 3t 

— 0.00006 cos 4t.-., 
to = - 0.12062 sin t -f 0.05079 sin 2t - 0.00184 sin 3t + 0.00096 sin 4t • • • . 

The infinitesimal body is in a line with the finite bodies and between them 
when T= 0. The value of r at this time is r(0)= 0.12667. The corre- 
sponding value given by Darwik is 0.1265, an agreement within the limits of 
accuracy of his computations. The infinitesimal body is in opposition at t ^tt, 
and we find for the radius at this time r(7r) = 0.11345. Darwin's value is 
0.1135. The agreements are sufficient to show that the analysis of this paper 
pertains to Darwin's "Satellites A." 

To get the coordinates in a retrograde orbit having the same sidereal period, 

it is necessary to change the sign of n in m = ?i'/(?i — 7i). In the example 

iust treated it is found that 

n' 61.383 



n ~ 421,383' 
whence m In the retrograde orbit having the same sidereal period is 
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n 
n 



m ■ 



1 + ^^ 
n 



61.383 - 1 

482.766 "~ 7.8648 ' 



The value of aja' is the same as before and the series give for the coordinates 

r= + 0.12412 - 0.00057 cos t - 0.00158 cos 2t 

- 0,00005 cos 3t - 0.00001 cos 4t + • • -, 
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w=+ 0.00820 sin t + 0.01654 sin 2t + 0.00047 sin 3t + 0.00010 sin 4t + • • • . 

These two orbits and the circular orbit of which they are the analytic continua- 
tion, together with the positions at intervals of j-V of the period, are given in 
the figure. 

Darwin found periodic planetary orbits corresponding to the " Satellites A." 
The analysis above gives these orbits within the limits of the convergence of the 
series if we put ^= 10, ^= I, M= 0, a' = 1, t? = S/{£J+S') = Jj . For 
the comparison with an actual example we shall take " Planet A " for the 
Jacobian constant equal to 40.0 (loc. cit., p. 225). The period given by Dar- 
win's computation is nT= 154° 13', whence 

"^ = ^t'ir = 2Mu = ^■^^^^^' ^' = ^•^^'^^^' -^^' = ^•^^^^- 

With these values of the parameters the series for the radius becomes 

r = 0.43373 + 0.00776 cos t - 0.01286 cos 2t - 0.00104 cos 3t • • .. 

When the infinitesimal body is between JiJ and S, and when it is in opposition 
to S, the values of r are respectively r = 0.42759 and r = 0.41415. Dar- 
win found in these respective cases r = 0.423 and r = 0.4140. These results 
agree within the limits of error of the computations and show that the " Planets 
A " belong to the type of orbits covered by the analysis of this paper. 

Darwin's " Satellites B" and "Satellites C" do not exist until the Jacobian 
constant has diminished nearly to 39.5. At first appearance they are coinci- 
dent, but with decreasing values of the Jacobian constant they gradually diverge. 
The computations of Darwin show that they are not isolated orbits existing 
only for special values of the parameter m. Suppose one of these orbits is 
known for m = m^. Unless the orbits "i?" and "C" coincide for this value 
of wi, the coordinates of an oi-bit near the known one may be expressed as power 
series in m — m^= X; if the orbits "5" and "C" coincide for this value of 
m, the coordinates of the two neighboring orbits may be expressed as power 
series in \'X and — Vx respectively. When X has one sign their coordinates are 
real, and when it has the other, they are conjugate complex quantities. Sup- 
pose the series converge so long as |X| <X|j. The coordinates of the orbits 
near the ones belonging to m = m^ + X may be expanded as power series in 
m —in^ — Xji = X' . If the complex orbits are in question there will be two, one 
with X and one with — X, whose coordinates will still be complex and conju- 
gate. They have the property of reentering after one synodic revolution. 

Suppose the two orbits which branch at m = wi„ are followed until 
m^ -f Xjj + ■ ■ • = . Their coordinates will remain conjugate through the whole 
series of steps, and the property of reentering after one synodic revolution will 
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be preserved. But for this value of m the disturbance vanishes and the differ- 
ential equations reduce to those which belong to the two body problem. The 
period remains real in t as well as in t. Consequently these orbits show that 
there are at least two conjugate complex solutions of the problem of two bodies 
whose period in t is real. 

Let us consider the differential equations expressed in t. From (6) we have 



n [ _ n \ 

m \ n ) 



and the differential equations may be expressed in terms of w'/w = M- The par- 
ameter enters in two ways, (a) independently of the trigonometric functions, and 
(6) under the trigonometric functions. We shall generalize the differential 
equations by distinguishing the it, according as it enters in one way or the other. 
Where it enters under the trigonometric functions it will be regarded as an 
assigned number which admits of no change in the development of the solutions ; 
where it enters in the other way it will be regarded as a parameter capable of 
taking all values sufficiently small, and it will be made to play the same role as 
m has in the preceding work. The solutions exist for all values of /a sufficiently 
small, but they belong to the physical problem only when the parameter /-i is set 
equal to the /* which occurs under the trigonometrical functions. The series of 
orbits is obtained by assigning a series of values to the /w. in the trigonometric 
functions. 

In the orbits in question the angular motion of the finite bodies is real; 
therefore n is real. Also the period of the infinitesimal body is real ; therefore 
/I and n are real. The orbit of the finite bodies is real ; therefore a is real 
and equal to unity with Darwin's choice of the linear unit. The parameter a 
is subjected only to the relation (7) 



,/ m \i ( E \i 
a ^ a I :i I I "^ j; I ( M being InflDltesimal ) . 



There are precisely three values of a satisfying this relation, one of which is real 
while the other two are conjugate complex quantities. All three satisfy the 
conditions of the two body problem equally. Likewise the periodic solutions of 
the three body problem exist for all three roots for all values of /i sufficiently 
small. Moreover these are the only solutions of the type under consideration 
in which the periods are real. In the three cases the coefficients of all powers 
of a I a' are precisely the same, and since the three values of aja' have the same 
modulus, the series have the same radius of convergence. Since these are the 
only solutions of the type in question, it follows that Darwin's orbits "Satel- 
lites B" and "Satellites C" are the analytic continuation of these complex 
orbits, and that his search for the satellite orbits was exhaustive. His work 
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shows the interesting and important result that the two imaginary orbits unite 
and become real (if the distribution of masses is 10 and 1) when m is a little 
less than I. 

The remaining problem respecting the series which have been developed is 
their true radius of convergence for a given distribution of mass between the 
finite bodies. For simplicity suppose the coordinates and components of velocity 
of the infinitesimal body are x^, •• ■, x^. Then the differential equations may 
be written 

dx, ^ , . 

'^=fA^i' •••,a'4> /^; (j=l,.--,4). 

The initial values of the x. are a, e, 0, 0. Then the solutions may be written 

(56) a!, = p^.(a, e, /i; «), 

which converge for any preassigned interval in t if the moduli of a, e and m are 
sufficiently small. Happel* has shown that they converge for the values of 
these parameters which belong to the moon's variational orbit. But this by no 
means proves the convergence of the series in fi which belong to the variational 
orbit. The parameters a and e are determined in terms of fi by the periodicity 
conditions (on two of the equations for symmetrical orbits) 

(57) ^.f a, e, /^; y j = 0, ^7 a, e, /i; y j == 0, 

where T is the period. These series are to be solved for a and e as power 
series in /t and the results substituted in (56). In general, new singularities 
are introduced by the conditions (67) and are carried back into (56). Conse- 
quently in discussing the convergency of the series for the variational orbit 
both the series (56) and the singularities introduced by (57) must be considered. 

Consider equations (57). For that value of /x for which two orbits coincide 
these equations have a double solution, and this value of /x is a branch point of 
their solutions. Therefore the periodic solutions in /x do not con verge, for 
values of /a whose moduli are greater than the modulus of that value of /u. for 
which two of the orbits coincide. Applying this to Darwin's computations, 
we find that the series for the complex orbits "5" and "C" do not converge 
for values of ft. beyond which these orbits unite and become real. The series 
for the orbit "J. " have the same radius of convergence as those for the complex 
orbits. Therefore the series for "J. " do not converge for values of /i greater 
than that for which " B" and "C" become real. Expressed in m, the series 
do not converge for | m| > J . In the orbit computed m = 1/5.86 , and an ex- 
amination of the coefficients shows that this is near the limit of practical useful- 

* Inaugural Dissertation, Gottlngen, 1900. 
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ness of the series. Many of Darwin's orbits belong to values of m for which 
the series diverge. 

Unfortunately It cannot be positively asserted that the series converge until 
the orbits " B " and " C " become real. In the first place the parameter ix must 
be considered as complex except as it appears under the trigonometric functions, 
where it must remain real. (The radius of convergence is entirely independent 
of the value of the parameter under the sine and cosine functions.) Now it 
might happen that there would be a branching of orbits for a smaller modulus 
of II in this complex realm. Such computations as those made by Darwin 
would fail to show it, for the orbits, real or imaginary, belong to the physical 
problem only when /u , as the argument of the series, has the same value as occurs 
under the trigonometric functions. 

The functions may branch in another way. An orbit which in general 
reenters only after two synodic revolutions may for a special value of ji coincide 
with one which reenters after one synodic revolution. At this point the orbit 
under consideration changes its stability.* Darwin's " Satellite A " loses its 
stability for m about 1/^.1. However, this seems not to have been because of 
a branching of the character in question, but because of a union of the orbit 
with the orbit of an " Oscillating Satellite a ." * But there is no earlier change 
of stability, at least for real values of /j, . 

The series do not have a fixed radius of convergence as t goes through a period. 
The convergence in question is a uniform convergence with respect to t over the 
period of motion, and in this sense we have spoken of the series as having radii 
of convergence. 

The problem of chief practical interest is the Lunar Theory. The distribu- 
tion of mass between the finite bodies assumed by Darwin, as compared with 
that occurring in the Lunar Theory, has decreased rj a little and has increased 
a/a for a given value of m about 30 fold. It seems probable that in the Lunar 
Theory the series converge for larger values of m than in Darwin's problem, 
and the computation seems to show that in the latter problem they converge for 
m = ^ . Although it can scarcely be doubted that the series converge for the 
values of the parameters which the Lunar Theory requires, a proof of their con- 
vergence for these values is yet wanting. However, we do know that there is a 
limit within which the series are perfectly valid. 
The University op Chicago. 

* Poinoaee's Les Methodes Nouvelles, vol. 3, art. 380. 



